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O ■ Abstract 



We predict geoneutrino fluxes in a reference model based on a detailed description of Earth's 



Qh. crust and mantle and using the best available information on the abundances of uranium, thorium, 



and potassium inside Earth's layers. We estimate the uncertainties of fluxes corresponding to 
the uncertainties of the element abundances. In addition to distance integrated fluxes, we also 
provide the differential fluxes as a function of distance from several sites of experimental interest. 
Event yields at several locations are estimated and their dependence on the neutrino oscillation 
parameters is discussed. At Kamioka we predict A^(U + Th) = 35 it 6 events for 10^^ proton yr and 
100% efficiency assuming s\t?{29) = 0.863 and 5m? = 7.3 x 10^^ eV^. The maximal prediction is 
55 events, obtained in a model with fully radiogenic production of the terrestrial heat flow. 
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I. INTRODUCTION 

By looking at antineutrinos from reactors, KamLAND jl| has confirmed the oscillation 
phenomenon previously discovered by SNO with solar neutrinos and has provided crucial 
information on the oscillation parameters. Putting together the results of solar and terres- 
trial experiments, the best fit is obtained at (5m^ = 7.3 x 10^^ eV^ and sin^(2^) = 0.863 j^. 
Since we know their fate from production to detection, neutrinos can now be used as physical 
probes. 

Furthermore, the detector is so pure and the sensitivity is so high that KamLAND will 
be capable of studying geoenutrinos, the antineutrinos originating from Earth's natural 
radioactivity. Indeed, from a fit to the experimental data the KamLAND Collaboration 
reported four events associated with ^ss-q ^^^^ with ^'^^Th decay chains. This result, 
obtained from an exposure of just 162 ton yr, provides the first insight into the radiogenic 
component of the terrestrial heat. KamLAND has thus opened a new window for studying 
Earth's interior and one expects more precise results in the near future from KamLAND 
and other detectors which are presently in preparation. 

The argument of geoneutrinos was introduced by Eder 0] in the 1960's and it was ex- 
tensively reviewed by Krauss et al. [5] in the 1980's. Raghavan et al. |^ and Rothschild et 
al. vA ^ remarked on the potential of KamLAND and Borexino for geoneutrino observations. 

n n n 

Fiorentini et al. j8|, |9|, [101] discussed the relevance of geoneutrinos for determining the radio- 
genic contribution to Earth's heat fiow and their potential for improving our knowledge of 
oscillation parameters, see also Ref. 

In preparation to the data which will become available in the near future, we present a 
systematic discussion of geoneutrinos, which incorporates the best geological and geochem- 
ical information on their sources and outlines the main uncertainties, so as to understand 
what can be gained from the study of geoneutrinos concerning both Earth's interior and 
neutrino properties. With this spirit, we shall consider the following points. 

(i) We provide a reference model that incorporates the best available knowledge for the 
distribution of U, Th, and K in Earth's interior. 

(ii) Within this model we predict neutrino fiuxes and signals for detectors at different 



^ We shall always refer to the version available as arXiv:nucl-ex/9710001 
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positions on Earth. 

(iii) We estimate uncertainties of neutrino fluxes and signals corresponding to uncertain- 
ties of the U, Th, and K distributions. 

II. THE REFERENCE MODEL: ELEMENT DISTRIBUTIONS AND GEONEU- 
TRINO FLUXES 

A global look at Earth's interior is useful before entering a detailed discussion on the 
element distributions. The amount of information which we (assume to) have on Earth's 
interior is somehow surprising, if one considers that the deepest hole which has ever been 
dug is probably only ten kilometers deep, a mere dent in planetary terms. 

The outer layer is the relatively thin crust, divided in two types, continental crust (CC) 
and oceanic crust (OC). The former averages 38 km in thickness, varying around the globe 
from 20 to 70 km, and it is made primarily of light elements such as potassium, sodium, 
silicon, calcium, and aluminium silicates. The oceanic crust is much thinner, from about 6 
to 8 km. 

Inside this crustal skin is Earth's mantle which is 2900 km deep overall. Largely made up 
of iron and magnesium silicates, the mantle as a whole accounts for about 68% of Earth's 
mass. One distingushes the upper mantle ^ (UM) from the lower mantle (LM), however, 
the seismical discontinuities between the two parts do not necessarily divide the mantle into 
layers. The main questions about the mantle are does it move as a single layer or as multiple 
layers? Is it homogeneous in composition or heterogeneous? How does it convect? These 
questions sound simple, but the answers are complex, possibly leading to more questions, 
see Ref. Q- 

Inside the mantle is Earth's core, which accounts for about 32% of Earth's mass. Based 
on comparison with the behavior of iron at high pressures and temperatures in laboratory 
experiments, on the seismic properties of the core, and on the fact that iron is the only 
sufficiently abundant heavy element in the universe, the core is generally believed to be 
made primarily of iron with small amounts of nickel and other elements. Over thirty years 
ago, however, it was suggested that a significant amount of potassium could be hidden in 

^ We shall define the upper mantle as the shallow mantle plus the transition region, i.e., the region below 
the crust down to 677 km 
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Earth's core, thus providing a large fraction of the terrestrial heat flow through K decay. 
This controversial possibility has been revived recently in Ref. jlS^ . 

Uranium, thorium, and potassium are lithophile elements, which accumulate in the con- 
tinental crust. Their abundance in the mantle is much smaller, however, the total amounts 
are comparable with those in the crust, due to the much larger mantle mass. The core is 
generally believed to contain negligible amounts of these elements. 

A global description of the present crust-plus-mantle system is provided by the bulk 
silicate earth (BSE) model, a reconstruction of the primordial mantle of Earth, subsequent 
to the core separation and prior to crust differentiation, based on geochemical arguments. 
In the BSE model the uranium abundance ^ is obseIU) = 2 x 10^*^, and one has Th/U 
= a(Th)/a(U) = 3.9 and K/U = a(K)/a(U) = 1.14 x 10^, where the quoted values are 
averages between different estimates, all consistent with each other to the level of 10% or 
better, see Table HI In the BSE model the total masses of uranium, thorium and potassium 
are thus M(U) = 0.81 x 10^^ kg, M(Th) = 3.16 x 10^^ kg, and M(K) = 0.49 x 10^^ kg. 

The equation relating masses and heat production is 

= 9.5M(U) + 2.7M(Th) + 3.6M(K) , (1) 

where H is in TW, M(U) and M(Th) are in units of 10^^ kg, and M(K) in units of 10^^ kg. 

In the BSE model, the contributed heat production rates are H{U) = 7.6 TW, if(Th) = 
8.5 TW, and -f^(K) = 1.8 TW, for a total of about one half of the observed terrestrial heat 
flow {He^^O TW). 



A. Uranium, thorium, and potassium distributions 

Our aim is to build a reference model (labeled as "ref"), which incorporates the best 
available knowledge of U, Th and K distributions inside Earth. Concerning Earth's crust, 
we distinguish oceans and seawater, the continental crust, subdivided into three sublayers 
(upper, middle, and lower), sediments and oceanic crust. All these layers have been mapped 
in Ref. , which provides values of density and depth over the globe on a grid with 2° steps. 
We distinguish next the upper mantle (extending down to about 600 km), the lower mantle 

^ We shall always refer to element abundances in mass and we remind the reader that the natural isotopic 
composition is ^ssy/U = 0.993, 232Th/Th = 1 and "Ok/K = 1.2 x IQ-*. 
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(down to about 2900 km), and the core, and use the prehminary reference earth model 
(PREM) for the values of the density at each depth, assuming spherical symmetry. 

For each component, one has to adopt a value for the abundances a(U), a(Th), and a(K). 
In the literature of the last twenty years one can find many estimates of abundances for the 
various components of the crust (OC, upper CC, lower CC, . . . ), generally without an error 
value (see Tables ITTl IIIH and II V|) . two classical reviews being in Refs. ly, till and a most 



useful source being provided by the GERM Reservoir database ^ ^ 

For the upper mantle we are aware of several estimates by Jochum et al. llOl , White l20| , 

n n n 

O'Nions and McKenzie [21|, Hofmann |22|, and Zartman and Haines |2J|. In this respect 
data obtained from material emerged from unknown depths are assumed to be representative 
of the average composition down to about 600 km. 

For each (sub)layer of the crust and for the upper mantle, we adopt as reference value for 
the uranium abundance a''*^^(U) the average of the values reported in Tables ITTl III H and IIVI 
Concerning Th and K, we observe that the abundance ratios with respect to uranium are 
much more consistent among different authors than the corresponding absolute abundances. 
We shall thus take the average of ratios and from these construct the reference abundances 
for thorium and potassium: 

a-f(Th) = (Th/U) a'-^f(U) and a''\K) = (K/U) a'^^V) . (2) 

For the lower mantle, where no observational data are available, we resort to the BSE model, 
which — we recall — describes the present crust-plus-mantle system based on geochemical 
arguments. 

The mass of each element (X = U, Th, K) in the lower mantle Mi,u{X) is thus obtained 
by subtracting from the BSE estimate the mass calculated for the crust and upper mantle: 

Mlm(X) = Mbse(X)-Mcc(X)-Moc(X)-Mum(X) . (3) 

Reference abundances for the lower mantle are then obtained by dividing these values by its 
mass Mlm = 2.9 x 10^^ kg. According to geochemical arguments, negligible amounts of U, 
Th and K should be present in the core. 

The resulting choice of input values for the reference model is collected in Tables InJ IIIH 
and IIVI Concerning this reference model, we remark the following points. 

(i) The uranium mass in the crust Mc(U) = 0.35 x 10^^ kg is mainly concentrated in 



the continental part. The oceanic crust contributes as httle as 0.005 x 10^ kg, since its 
impoverished by a factor of 20 and it is much thinner than the continental crust. 

(ii) The estimated uranium mass in the upper mantle is about one sixth of that in the 
crust, whereas the lower mantle contains about as much uranium as the crust. 

(iii) Note that in this refererence model, constructed so as to satifisy the BSE constraint 
(121), mantle depletion (with respect to BSE) extends to the lower mantle. 

(iv) Similar considerations hold for thorium and potassium. 

B. The reference fluxes 

For each element X the produced ^ antineutrino fluxes at position r are defined as 



where r is the lifetime, is the atom mass, n is the number of antineutrinos per decay chain, 
the integral is over the volume of the earth, p is the local density, and ax is the abundance of 
the element X. We have evaluated the produced fluxes at several sites on the globe within 
the reference model (a = a'^'^^). We concentrate here on a few locations of specific interest, 
see Tables El EH and IVTTl 

(i) for the Kamioka mine, where the KamLAND detector is in operation, we predict an 
uranium flux $[/ = 3.7 x 10^ cm~^ s~^, a comparable flux from thorium and a fourfold 
flux from potassium. Within the reference model, about 3/4 of the flux is generated from 
material in the crust and the rest mainly from the lower mantle. 

(ii) At Gran Sasso laboratory, where Borexino i24^ is in preparation, we predict an ura- 
nium flux $u = 4.2 X 10^ cm~^ s~^, this larger flux arising from a bigger contribution of the 
surrounding continental crust. Thorium and potassium fluxes are correspondingly rescaled. 

(iii) At the top of Himalaya, a place chosen so that the crust contribution is maximal, we 
find the maximum uranium flux $u = 6.7 x 10^ cm~^ s^^. The crust contribution exceeds 



(iv) On the Hawaii, a site which minimizes the crust contribution, we find $u = 1.3 x 
10^ cm~^ s~^, originated mainly from the mantle. 

The produced fluxes are calculated ignoring oscillations, which will be discussed later. 




(4) 



90%. 
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These computed reference fluxes are generally larger than those of Rothschild [3], by a 
factor of order 30-50 %. This arises from several differences in the approaches. 

(i) We have used a more recent and detailed map of Earth's crust: the grid is denser and 
several layers are distinguished. 

(ii) We have a more detailed model for the mantle, correponding to the PREM density 
proflle. 

(iii) Most important, our reference values for the abundances in the continental crust are 



larger than that used in Ref. [7[. As an example, Rothschild et al. use for the continental 
crust acc(U) = 0.91 ppm from a classical review paper of 1985 16j. Our reference model, 
when averaged over the different sublayers, yields acc(U) = 1.5 ppm. This larger value arises 
from taking into account recent data, which are all higher than those quoted in Ref. Q]. 

The produced fluxes are computed ignoring the effect of oscillations, which depends on 
the distance R between the detector and the source. For taking into account this effect, and 
also in view of understanding which portion of Earth can be accessed with a geoneutrino 
detector, it is useful to introduce quantities which contain more detailed information. 

The differential fluxes f{R) are obtained by grouping together all the sources which lie 
at the same distance R from the detector 



Note that fx{R) actually depends also on the detector position r and just for semplicity of 
notation we drop this variable. 

The cumulated fluxes (f){R) are deflned as 



They represent the cumulative effect of all sources within a distance R from the detector: 
the total produced fluxes of Eq. (jH) are clearly $x = 0x(2-Re). 

These quantities are plotted in Figs. ^ and |21 for the four sites (we only show the uranium 
contribution, the shapes of the other contributions being similar). We remark that we 
have been using an average density approximation, which presumably breaks down near the 
detector, where one should resort to a detailed geological study of the surroundings. From 
Fig. 121 one sees that in our model the region within 30 km from Kamioka or Gran Sasso 
originates about 15% of the total produced flux. Concerning the region where most of the 






(6) 
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flux is generated, one sees again from Fig.|21tliat 50% of the produced flux originated witliin 
400 km (800 km) from Kamioka (Gran Sasso). 

In Tables rVllll IIXI and 1X1 we present the numerical values of the contribution to fx{R) 
from the crust at Kamioka and Gran Sasso and that from the mantle (the assumed spherical 
symmetry of the mantle implies the same contribution at any site). These data will be useful 
for a detailed analysis of future experiments devoted to the study of geoneutrinos, in order 
to take into account the distance dependence of the survival probability. 



III. THE UNCERTAINTIES OF THE REFERENCE MODEL 

The fluxes of the reference model correspond to the best available knowledge about the 
crust and the interior of Earth, as derived from observational data and geochemical infor- 
mation on the global properties. An estimate of the uncertainties of the predicted fluxes is 
clearly useful. 

Since the abundance ratios look relatively well determined, we concentrate on the un- 
certainties of the uranium abundances in the different layers and propagate them to the 
other elements. For the reference model, we have Mcc(U) = 0.345 x 10^'' kg, Moc(U) = 
0.005 X 10^'' kg, the total mass of CC being Mcc = 2.234 x 10^^ kg. According to our model, 
the average uranium abundance in the CC is thus acc(U) = 1.54 x 10~^. We determine 
lower and upper limits by observing that the range of estimated uranium abundances is 
between 0.91 x 10~ 6 Q and 1.8 X 10-6 Q 

low: acc(U) = 0.9 x 10"^; acc(Th) = 3.7 x 10"^; acc(K) = 0.94 x 10"^ , 
high: acc(U) = 1.8 x 10"^; acc(Th) = 7.6 x 10"^; acc(K) = 1.97 x 10"^ . 

We remark that there is an overall uncertainty of a factor 2 concerning the total amount of 
radioactive materials in the crust. 



For the upper mantle, we take as extrema the two values known to us 19|,|23| for uranium 
and we deduce thorium and potassium by rescaling 

low: auM(U) = 5 x 10"^; auM(Th) = 13 x 10"^; auM(K) = 6 x 10"^ , 
high: auM(U) = 8 x 10"^; auM(Th) = 21 x 10"^; auM(K) = 9.6 x 10"^ . 

Such a small uncertainty is perhaps optimistic, however, it is not influential for the future 
discussion in view of the relatively small amounts contained in the upper mantle. 
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We remind the reader that no observational information is available for the lower mantle. 
For building a minimal model, we assume that the mantle is fully mixed and use for the 
whole mantle the lowest values estimated from samples coming from the upper mantle. 

A maximal model can be obtained by assuming that the terrestrial heat is fully accounted 
by radiogenic production. This can be otained by keeping the BSE abundance ratios fixed 
and rescaling the total masses to M(U) = 1.67 x 10^'' kg, M(Th) = 6.5 x 10^'^ kg, and 
M(K) = 1.9 X 10^^ kg. ^ A natural implementation is obtained by choosing for the crust 
and upper mantle the highest observational estimates and placing the remaining mass in 
the lower mantle. 

All this leads to 

low: ai^uiU) = 5 x lO'^; aLAf(Th) = 13 x lO'^; aLAf(K) = 6 x 10"^ , 
high: aLAf(U) = 40 X IQ-^; aLAf(Th) = 156 x 10"^; a^Af (K) = 45.6 x 10^^ . 

The corresponding low and high estimates of the produced fluxes are also shown in Tables |3 
IVIl and I VI II for a few locations. 

In view of assigning an uncertainty to the fluxes of the reference model one can take two 
different approaches. 

(a) A conservative estimate: the error assigned to the reference value is half of the 
difference between the high and low estimates A$cons = ('^'high — *I'iow)/2. 

(b) A statistical estimate: one assumes that the full range of calculated fluxes represents 
a ±3cr interval. ® In this way one obtain a conventional la error A$ = ($high — *^'iow)/6. 

The relative uncertainties of the fluxes are reported in Table IXII They are the same 
(and fully correlated) for all elements, the l-cr error being about 15%, at Kamioka and Gran 
Sasso. At Hawaii, where the mantle contribution is dominant, the error is much larger, as a 
consequence of the large uncertainties of the lower mantle's composition. 

When using these errors, one has to remark that uncertainties associated with abundances 
in the crust and in the upper mantle are deduced from the spread of observational data, 
whereas the estimates for the lower mantle, which cannot be accessed by observations, 
completely rely on theoretical arguments. In addition, one should also take into account the 
detailed geological structure around the detector for more precise flux estimates. 

^ Clearly this model does not satisfy the BSE constraint on the total U, Th, and K masses in the Earth. 
^ If unhappy with this conventional assumption, the reader can rescale a. 
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IV. FROM FLUXES TO SIGNALS AND DETECTORS 



Geoneutrinos can be detected by means of inverse beta reactions 

z7e + (^,A)^e+ + (Z-l,A) , (7) 

where the positron kinetic energy T is related to the antineutrino energy E hj T = E — Eq, 
with Eq = mz-i + rrif, — rriz- The differential event yield as a function of T is given by 



dN /'2-Re 

— = Nzta{E)Y,wx{E) dRfx{R)Pee{E,R) , (8) 

X 



where Nzt is the exposure (number of target nuclei times the live time), cT(i?) is the cross 
section of reaction T = E — Eq and the integral is over the distance R from the detector. 
The survival probability of z/g produced at distance R with energy E is 

P^E, R) = l- sm'{2e) sin^ 1^^) ■ (9) 

For each element, the differential produced flux fx(R) is defined in Eq. (jSj), wx{E) is the 
energy spectrum of the Ue from the decay chain 26| of element X and normalized to 1, 
Jq°° dEwx{E) = 1. For simplicity we neglect the finite energy resolution of the detector and 
assume 100% detection efficiency. 

Another interesting observable is the total geoneutrino yield 

where T^ax is the maximal positron energy. 

The classical approach to low energy antineutrino detection is by using hydrogen com- 
pounds as target, by means of i^e + p — > e+ + n. Since Eq = rrin + nie — nip = 1.804 MeV, 
this reaction is suitable for antineutrinos from uranium and thorium progenies (-Emax = 3.26 
and 2.25 MeV, respectively), whereas antineutrinos from potassium (-Emax = 1-31 MeV) are 
below threshold. 



^ A frequently used variable is the visible energy E^is — T + 2™^ which is the energy released in the slowing 
down and subsequent annihilation of the positron. 
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A. Total yields 



We discuss first the total geoneutrino yield A^, which is experimentally more accessible 
than the differential spectrum. In view of the structure of the survival probability, see 
Eq. (0), it can be written as 

iV = iV,4l-sin2(2e)x(W)] , (11) 

where A^no is the yield for no oscillation. 

The function x depends on the uranium and thorium distributions inside Rarth and on 
the detector position. Obviously x tends to (1/2) for small (large) values of 5w?. We have 
computed x in the reference model for some sites of interest, see Fig. El At all locations and 
for 5m? > 4 X 10"^ eV^, the function x differs from its asymptotic value by less than 2%. 
Using the asymptotic value of the survival probability and the best fit value of the mixing 
angle jsj, one finds 

N = Nr,o [I - 0.5 sm^{29)] = 0.57Ar„o . (12) 

In Fig. 0] we show the relative contributions of different distances to the total yield: for the 
most interesting values of the region within 30 km from Kamioka contributes about 
15% of the total. The no oscillation yield A^no is determined in terms of the total produced 
fluxes from uranium and thorium 

N^o = 13.2$u + 4.0$Th (13) 

for an exposure of 10^^ proton yr with fluxes $ in units of 10^ cm^^ s^^. 

The no oscillation yields, calculated with the fluxes of the reference model, are shown 
in Table IXIII In the same table we also present the estimated la errors, obtained by 
propagating those on the produced fluxes (which are dominant over the other uncertainties 
from cross sections, decay spectrum, etc.) and the minimal and maximal predictions. 

For the Kamioka site the prediction of the reference model (normalized ^ to 10^^ proton yr 
and 100% efficiency) is A^no = 61 in good agreement with the "best model" of Refs. [sl llQ]. 



^no = 67 TNU, in between the values of Ref. [3], A^no = 43 TNU, and of model lb of 



^ It is useful to introduce a terrestrial neutrino unit (TNU) for event rates, defined as one event per 10'^^ 
target nuclei per year, or 3.17110"*° per target nucleus. This unit is analogous to the solar neutrino 
unit (SNU) 0. 
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Ref. p, A^no = 75 TNU. An experimental value for A^no = 156 TNU can be deduced from 
the nine geoevents reported by KamLAND, assuming Pee = 0.57. All the above predictions 
are consistent with the experimental result within its statistical error (about 60% [^). 

The total yields predicted in our reference model for a number of locations are presented 
in Fig. El We remind the reader that geoneutrino fluxes are superimposed to the low-energy 
tail of antineutrinos from nuclear reactors, which can provide in this respect an important 



background, as first pointed out by Lagage 



28l. This effect is clearly dependent on location 



and it has been extensively discussed in Ref. [lO|. In particular, the event yield from reactors 
has been estimated as about 300 TNU (no oscillations) at Kamioka and about 70 TNU at 
Gran Sasso. 



B. Event spectra 

A more detailed information is contained in the event spectrum dN/dT and a relevant 
question is whether the spectrum is deformed because of oscillations. From Eqs (jHl) and 
the event distribution with energy can be written as 

^=(^) [l-sin^(2^)^(T,W)] , (14) 

where T is the positron kinetic energy. 

The no-oscillation spectrum dN^^/ dT is shown in Fig. IHl for Kamikande. The function 
■ip{T, 5rn?) represents the modification to the event spectrum due to oscillations. It is plotted 
for Kamioka for a few values of (5m^ in Fig. [3 One sees that oscillations produce a moderate 
distortion for the two smallest values of 5m? and the distortion is negligible for the largest 
values of dm?. 



V. CONCLUDING REMARKS 



We summarize here the main points of this paper 

(i) We have provided a reference model for the produced fluxes of geoneutrinos, estimating 
its uncertainties in view of available data and geochemical inferences about U, Th, and K 
distribution in Earth's interior. When normalized to an exposure of 10'^^ proton yr, an 
averaged survival probability {Pee) = 0.57 and a 100% detection efficiency, the predicted 
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events for KamLAND are 



A^(U) = 28±4.7 , A^(Th) = 7±1.2 . (15) 

Errors have been estimated so as correspond to la confidence level and are (almost) com- 
pletely correlated: 

A^(U + Th) = 35 ± 6 . (16) 

(ii) Concerning the estimated errors, we remark that uncertainties associated with abun- 
dances in the crust and in the upper mantle are deduced from the spread of observational 
data, whereas the estimates for the lower mantle, which cannot be accessed by observations, 
completely rely on theoretical arguments. In addition, one should also take into account the 
detailed geological structure around the detector for more precise flux estimates. 

(iii) We have also investigated extreme models, corresponding the the minimal and max- 
imal amounts of U and Th which could be present on Earth. At KamLAND we predict 

^low^U + Th) = 29 and N^'^^{\J + Th) = 74 . (17) 

In these two extreme models U, Th, and K, in the BSE proportions, produce a radiogenic 
heat ifrad = 9 and 40 TW, respectively. If experimental results come out close to the 
minimum, then uranium and thorium provide a minor contribution to Earth's energetics: 
either Earth's heat flow is mainly non radiogenic or a significant amount of potassium has 
to be hidden in Earth's interior. If values near to the maximal are found from experiments, 
then radiogenic contribution is the main supply of Earth's heat flow, and one can exclude 
models where significant amounts of potassium are hidden in Earth's core. 

(iv) Predictions for detectors at several locations are also given, see Table IXTTl and Fig. 
We remark that a detector located far from the continental crust could provide significant 
information on the structure of the mantle, particularly when compared with data from 
detectors at sites where (as in KamLAND and Borexino) the contribution of Earth's crust 
is important. 
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TABLE I: Abundances in the bulk silicate Earth model. 



a(U) Th/U K/U Remarks 

2.1 X 10~^ 4.0 1.14 X 10^ [22] 

2.3 X 10-8 [211] 

2.0 X 10-^ 4.0 1.27 X 10^ [22] 

1.8 X 10-8 3.6 1.0 X 10^ [16] 

2.0 X 10-8 3.9 1.14 X 10^ average 



TABLE II: Uranium abundances in Earth's interior. 



Layer 




Available data 


Adopted value Remarks 








a(U) 








Oceans & Seawater 




3.2 X 10-9 


3.2 X 10-9 


m 




Sediments 




1.68 X 10-^ 


1.68 X 10-6 


m 




Upper CC 


(2.2 


;2.4 ;2.5 ;2.8) x IQ-^ 


2.5 X 10-6 


Average of [12], 


m 


Middle CC 




1.6 X 10-6 


1.6 X 10-6 


m 




Lower CC 


(0.20 ; 


0.28 ;0.93 ;1.1) x IQ-* 


' 0.62 X 10-6 


Average of [34], [16], [17], 


[25] 


Oceanic crust 




0.1 X 10-6 


0.1 X 10-6 


\16] 




Upper mantle 




(5 ;8) X 10-9 


6.5 X 10-9 


Average of [19], [23] 




Lower mantle 






13.2 X 10-9 


From Eq. (jH)) with 

aBSE(U) = 2 X 10-8 




Core 
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TABLE III: Thorium abundances in Earth's interior. 



Layer 


Available data 


Average 


Adopted value Remarks 






Th /TT 




(Th/U) 


Lt I -1 11 J 






Oceans & Seawater 


n 
u 







u 






Sediments 


A 1 1 




4.11 








Upper CC 


3.8 ;3.8 ;3.9 ; 


4.1 


3.9 


9.8 X lO"'^ 


Average of [34], [33], [33], 


[17] 


Middle CC 


3.8 




3.8 


6.1 X lO-'^ 


m 




Lower CC 


3.8 ; 6.0 ;7.0 ; 


7.1 


6 


3.7 X 10-6 


Average of [16], [34], [25], 


[17] 


Oceanic crust 


2.2 




2.2 


0.22 X 10^6 


[16] 




Upper mantle 


2.58 ;2.63 ;2.7 


;2.73 


2.66 


17.3 X 10^9 


Average of [20], [21j, [22], 


[23] 


Lower mantle 








52.0 X 10-9 


From Eq. © with 





aBSE(Th) = 7.8 X 10-8 

Core 



TABLE IV: Potassium abundances in Earth's interior. 



Layer 



Available data 



Average Adopted value Remarks 





(K/U) X 10-4 


(K/U) X 10~4 




Oceans & Seawater 


12.5 


12.5 


4.0 X lO-'' 


Sediments 


1.0 


1.0 


1.7 X 10-2 


Upper CC 


0.99 ;1.0 ;1.03 ; 1.10 


1.03 


2.57 X 10-2 


Middle CC 


1.04 


1.04 


1.67 X 10-2 


Lower CC 


1; 1.2; 1.4 


1.2 


0.72 X 10-2 


Oceanic crust 


1.25 


1.25 


0.125 X 10-2 


Upper mantle 






0.78 X 10-4 


Lower mantle 






1.6 X 10-4 



m 

[32] 

Average of [16], [17], [33], [33] 

m 

Average of [Ij, [2^, [VJj 
[16] 



aBSE(K) = 2.32 X 10" 



Core 
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TABLE V: Uranium: masses, radiogenic heat, and predicted fluxes. Units are 10^^ kg, TW and 
10^ cm~^ s~^, respectively. The reference values, lower and upper limits are labeled as ref, low, 
and high, respectively. Crust summarizes CC and OC; UM (LM) denotes upper (lower) mantle. 









l~rilTlRlflVR 

-■- j-x-LiJ.cvJ.CA' y cti 


Ciran Sasso 


Kamioka 










33° N 85° E 


42° N 14° E 


36° N 137° E 


20° N 156° W 




IVl \\J J 


n yyj) 










r^m^t; low 


0.206 


1.960 


3.337 


1.913 


1.594 


0.218 


V^'X LlOU J- CJ. 


0.353 


3.354 


5.710 


3.273 


2.727 


0.373 


liTiicj'l" ni o" n 


413 


3 920 


6 674 


3 826 


3 187 


436 


TIM low 






146 


146 


n 146 


146 


UM ref 


Ofi2 


.591 


1 89 


189 


1 89 


189 


UM high 


0.077 


0.727 


0.233 


0.233 


0.233 


0.233 


LM low 


0.147 


1.399 


0.288 


0.288 


0.288 


0.288 


LM ref 


0.389 


3.695 


0.760 


0.760 


0.760 


0.760 


LM high 


1.177 


11.182 


2.299 


2.299 


2.299 


2.299 


Total low 


0.401 


3.814 


3.770 


2.346 


2.027 


0.651 


Total ref 


0.804 


7.639 


6.659 


4.222 


3.676 


1.322 


Total high 


1.666 


15.828 


9.206 


6.358 


5.720 


2.968 
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TABLE VI: Thorium: masses, radiogenic heat, and predicted fluxes. Units are 10^^ kg, TW and 
10^ cm~^ s~^, respectively. The reference values, lower and upper limits are labeled as ref, low, 
and high, respectively. Crust summarizes CC and OC; UM (LM) denotes upper (lower) mantle. 









Hima.l3.ycL 


Ciran Sasso 


Ka.miokB. 


HcLWciii 








33° N 85° E 


42° N 14° E 


36° N 137° E 


20° N 156° W 




1V± I J. 11 I 












r^m^t; low 


0.838 


2.263 


2.972 


1.714 


1.420 


0.180 


V^'X LlOU J- CJ. 


1.450 


3.915 


5.141 


2.964 


2.456 


0.311 


liTiicj'l" ni o" n 


1.722 


4 64Q 


6 105 


3 520 


2 Q16 


370 


TIM low 


0.124 




VJ .VJO(J 


083 


083 


083 


UM ref 


0.166 


0.447 


0.111 


0.111 


0.111 


0.111 


UM high 


0.207 


0.558 


0.138 


0.138 


0.138 


0.138 


LM low 


0.383 


1.034 


0.165 


0.165 


0.165 


0.165 


LM ref 


1.532 


4.135 


0.658 


0.658 


0.658 


0.65 


LM high 


4.590 


12.393 


1.973 


1.973 


1.973 


1.973 


Total low 


1.346 


3.633 


3.220 


1.961 


1.668 


0.428 


Total ref 


3.147 


8.497 


5.910 


3.733 


3.225 


1.080 


Total high 


6.519 


17.600 


8.216 


5.631 


5.028 


2.481 
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TABLE VII: Potassium: masses, radiogenic heat, and predicted fluxes. Units are 10 kg, TW, 
and 10^ cm~^ s~^, respectively. The reference values, lower and upper limits are labelled as ref, low, 
and high, respectively. Crust summarizes CC and OC; UM (LM) denotes upper (lower) mantle. 









l~rilTlPl,lflVR 

J. j.xj.j.j.cvJ.CA' y cij 


Oran Sasso 


Kamioka 


HcLWciii 








33° N 85° E 


42° N 14° E 


36° N 137° E 


20° N 156° W 
















r^Tn<3'i' low 


f) 21f) 


n 757 


12.429 


7 126 


5 Q41 


851 


V^'X LlOU J- CJ. 


0.367 


1.321 


21.684 


12.432 


10.366 


1.485 


liTiicj'l" ni o" n 


0.441 


1 587 


2fi 048 


14 934 


12.451 


1 784 


TIM low 








634 


n 634 


634 


UM ref 


075 




824 


824 


n 824 


824 


UM high 


0.092 


0.331 


1.015 


1.015 


1.015 


1.015 


LM low 


0.177 


0.636 


1.254 


1.254 


1.254 


1.25 


LM ref 


0.471 


1.697 


3.343 


3.343 


3.343 


3.34 


LM high 


1.344 


4.838 


9.534 


9.534 


9.534 


9.534 


Total low 


0.444 


1.600 


14.317 


9.014 


7.829 


2.739 


Total ref 


0.913 


3.287 


25.852 


16.600 


14.533 


5.652 


Total high 


1.877 


6.756 


36.596 


25.482 


23.000 


12.332 
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TABLE VIII: Differential produced fluxes: the contributions from the crust at Kamioka. The 



distance R is in km, fx in cm ^ s 



R 


/u 




/k 


10 


1.86E-01 


1.61E-01 


6.94E-01 


20 


2.29E-01 


1.97E-01 


8.55E-01 


30 


2.01E-01 


1.75E-01 


7.57E-01 


40 


1.59E-01 


1.45E-01 


6.07E-01 


50 


1.23E-01 


1.12E-01 


4.68E-01 


60 


9.86E-02 


9.04E-02 


3.76E-01 


70 


8.34E-02 


7.65E-02 


3.18E-01 


80 


7.51E-02 


6.87E-02 


2.86E-01 


90 


6.62E-02 


6.06E-02 


2.52E-01 


100 


5.57E-02 


5.11E-02 


2.12E-01 


200 


2.31E-02 


2.12E-02 


8.82E-02 


300 


8.15E-03 


7.39E-03 


3.12E-02 


400 


5.24E-03 


4.74E-03 


2.01E-02 


500 


3.68E-03 


3.31E-03 


1.41E-02 


600 


2.61E-03 


2.35E-03 


l.OOE-02 


700 


2.47E-03 


2.23E-03 


9.50E-03 


800 


2.53E-03 


2.29E-03 


9.68E-03 


900 


2.94E-03 


2.67E-03 


1.13E-02 


1000 


2.88E-03 


2.61E-03 


l.lOE-02 


2000 


1.32E-03 


1.20E-03 


5.06E-03 


3000 


1.08E-03 


9.72E-04 


4.11E-03 


4000 


1.05E-03 


9.51E-04 


4.01E-03 


5000 


7.44E-04 


6.75E-04 


2.84E-03 


6000 


4.88E-04 


4.40E-04 


1.86E-03 


7000 


4.28E-04 


3.86E-04 


1.64E-03 


8000 


2.99E-04 


2.69E-04 


1.14E-03 


9000 


2.53E-04 


2.27E-04 


9.67E-04 


10000 


2.19E-04 


1.98E-04 


8.41E-04 


11000 


2.16E-04 


1.96E-04 


8.28E-04 


12000 


1.40E-04 


1.24E-04 


5.35E-04 
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TABLE IX: Differential produced fluxes: the contributions from the crust at Gran Sasso. The 



distance R is in km, fx in cm ^ s 



R 


/u 


f'Yh 


/k 


10 


1.48E-01 


1.29E-01 


5.48E-01 


20 


2.11E-01 


1.82E-01 


7.86E-01 


30 


1.80E-01 


1.59E-01 


6.79E-01 


40 


1.40E-01 


1.28E-01 


5.34E-01 


50 


1.12E-01 


1.03E-01 


4.26E-01 


60 


8.94E-02 


8.21E-02 


3.41E-01 


70 


7.66E-02 


7.04E-02 


2.92E-01 


80 


6.59E-02 


6.05E-02 


2.51E-01 


90 


5.92E-02 


5.43E-02 


2.25E-01 


100 


5.22E-02 


4.79E-02 


1.99E-01 


200 


2.30E-02 


2.11E-02 


8.75E-02 


300 


1.31E-02 


1.20E-02 


5.02E-02 


400 


1.14E-02 


1.04E-02 


4.34E-02 


500 


9.83E-03 


8.95E-03 


3.74E-02 


600 


7.52E-03 


6.81E-03 


2.86E-02 


700 


5.98E-03 


5.43E-03 


2.27E-02 


800 


5.01E-03 


4.56E-03 


1.91E-02 


900 


4.95E-03 


4.52E-03 


1.88E-02 


1000 


5.12E-03 


4.68E-03 


1.95E-02 


2000 


2.98E-03 


2.71E-03 


1.13E-02 


3000 


1.60E-03 


1.45E-03 


6.08E-03 


4000 


1.22E-03 


l.llE-03 


4.66E-03 


5000 


7.65E-04 


6.91E-04 


2.91E-03 


6000 


5.98E-04 


5.42E-04 


2.28E-03 


7000 


5.66E-04 


5.14E-04 


2.16E-03 


8000 


4.44E-04 


4.02E-04 


1.69E-03 


9000 


2.20E-04 


1.97E-04 


8.41E-04 


10000 


8.20E-05 


7.18E-05 


3.19E-04 


11000 


1.61E-04 


1.46E-04 


6.20E-04 


12000 


1.27E-04 


1.14E-04 


4.88E-04 
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TABLE X: Differential produced fluxes: the contributions from the mantle. The distance R is in 



km, fx in cm ^ s 



R 


fu 




/Th 




/k 


10 


O.OOE+00 




O.OOE+00 




O.OOE+00 


20 


O.OOE+00 




O.OOE+00 




O.OOE+00 


30 


O.OOE+00 




O.OOE+00 




O.OOE+00 


40 


O.OOE+00 




O.OOE+00 




O.OOE+00 


50 


1.62E-04 




9.48E-05 




7.05E-04 


60 


2.91E-04 




1.70E-04 




1.27E-03 


70 


3.77E-04 




2.21E-04 




1.64E-03 


80 


4.38E-04 




2.57E-04 




1.91E-03 


90 


4.84E-04 




2.83E-04 




2.11E-03 


100 


5.19E-04 




3.04E-04 




2.26E-03 


200 


6.64E-04 




3.89E-04 




2.89E-03 


300 


7.08E-04 




4.14E-04 




3.08E-03 


400 


7.31E-04 




4.28E-04 




3.19E-03 


500 


7.53E-04 




4.41E-04 




3.28E-03 


600 


7.71E-04 




4.51E-04 




3.36E-03 


700 


8.49E-04 




5.32E-04 




3.70E-03 


800 


9.83E-04 




6.74E-04 




4.30E-03 


900 


1.09E-03 




7.83E-04 




4.75E-03 


1000 


1.17E-03 




8.70E-04 




5.12E-03 


2000 


1.49E-03 




1.22E-03 




6.56E-03 


3000 


1.48E-03 




1.24E-03 




6.51E-03 


4000 


l.llE-03 




9.27E-04 




4.88E-03 


5000 


8.88E-04 




7.41E-04 




3.90E-03 


6000 


7.40E-04 




6.17E-04 




3.2r5E-03 


7000 


6.34E-04 




5.29E-04 




2.79E-03 


8000 


5.54E-04 




4.63E-04 




2.44E-03 


9000 


4.93E-04 




4.11E-04 




2.17E-03 


10000 


4.24E-04 




3.53E-04 




1.86E-03 


11000 


2.35E-04 




1.91E-04 




1.03E-03 


12000 


6.10E-05 




4.11E-05 




2.66E-04 


TABLE XL Fractional uncertainties of the produced fluxes. 


A$/$ (%) 


Himalaya 


Gran Sasso 


Kamioka 


Hawaii 


conventional 1 a 


14 


16 




17 


29 
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TABLE XII: Total yields. A^no is the total number of gcocvcnts (U + Th) in the absence of 
oscillations predicted from the reference model for 10^^ proton yr (or in TNU) and AA^no is the 
"Icr" error. (A^o^^) is the minimal (maximal) prediction. For SmP' > 4 x 10~^ eV^ the 

gcoovc-iit yi(^ld is A' = A^o [l - J) nin^ (20)] . 



Location 


Nno 


AAr„o 


Ajlow 

no 




Baksan 


91 


13 


51 


131 


Hawaii 


22 


6 


10 


49 


Himalaya 


112 


15 


63 


154 


Homestake 


91 


13 


51 


130 


Kamioka 


61 


10 


33 


96 


La Palma 


37 


8 


19 


67 


LGS 


71 


11 


39 


106 


Pyhasalmi 


92 


13 


51 


131 


Sudbury 


87 


13 


48 


125 


Yucca Mountain 


70 


11 


38 


106 
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FIG. 1: Differential produced flux from uranium as a function of the distance R from the detector. 
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FIG. 2: Cumulated produced flux from uranium as a function of the distance R from the detector. 
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FIG. 3: Dependence of the yield on on 6m^. The figure shows the function x = (-^no — 
N)/ [(iVno sin^(26'))] , see Eq. ((TT|) . for four locations with dm^ in units of 10^^ eV^. Solid (dashed, 
dotted, dot-dashed) line applies to Kamioka (LNGS, Hawaii, Himalaya). 
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FIG. 4: Contributed signal as a function of distance. The percentage contribution to the event 
yield at Kamioka originating from sources within R is shown for the indicated values of dm^ in 
units of 10^^ cV^ at fixed sin^(20) = 0.863. The percentage contributed neutrino flux without 
oscillation is also shown for comparison. 
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Q I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I II 

1 1.25 1.5 1.75 2 2.25 2.5 



Evis [MeV] 

FIG. 6: Event spectrum as function of the visible energy i?vis = T + 2me in MeV. The spectrum is 
calculated for the U/Th flux ratio expected at Kamioka with no oscillation and the normalization 
is arbitrary. 
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FIG. 7: Spectrum deformation. The function defined in Eq. (|14|) . as function of the visible 
energy E'vis = T + 2me in MeV for four values of dm?: 1 x 10^^ eV^ (dash hue), 3 x 10^^ eV^ 
(dot-dash line), 7.3 x 10"^ eV^ (sohd hue), and 20 x 10"^ eV^ (dot hue). 
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